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Finite temperature Casimir theory of the Dirichlet scalar field is 
developed, assuming that there is a conventional Casimir setup in 
physical space with two infinitely large plates separated by a gap R 
and in addition an arbitrary number q of extra compacified dimen- 
sions. As a generalization of earlier theory, we assume in the first part 
of the paper that there is a scalar 'refractive index' N filling the whole 
of the physical space region. After presenting general expressions for 
free energy and Casimir forces we focus on the low temperature case, 
as this is of main physical interest both for force measurements and 
also for issues related to entropy and the Nernst theorem. Thereafter, 
in the second part we analyze dispersive properties, assuming for sim- 
plicity q = 1, by taking into account dispersion associated with the 
first Matsubara frequency only. The medium-induced contribution to 
the free energy, and pressure, is calculated at low temperatures. 

PACS numbers: ll.10.Kk, ll.10.Wx, 42.50.Lc 

1 Introduction 

Consider two infinite parallel plates separated by a gap R, and assume that 
the field between the plates, as well as on the outside, is a scalar field obeying 
Dirichlet boundary conditions at z = and z = R (z is the direction normal 
to the plates). In order to get a Casimir configuration we have to take into 
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account the field outside the plates also. Let one of the plates, the right one 
at z = R, say, be denoted a "piston". Now generalize the situation such 
that p spatial dimensions are envisaged, together with q extra compactified 
dimensions. The spacetime dimension is thus D — p + q + 1. We are led to 
a Casimir piston model in which spacetime is flat. This model has attracted 
considerable attention in the recent literature PHEJUSISHEJEJCEnj 
HH H2l H3l 1131 HSl IIEl IIZl IIHl dSl I201 IS]. One reason for the 

current interest is obviously the mathematical elegance of the formalism; the 
efficiency of regularization procedures like the zeta function regularization is 
quite striking. Typical for this kind of theories is that the field energy can be 
expressed in terms of Epstein-like zeta functions. The Casimir force between 
the two plates in physical space follows by taking the derivative of the energy 
(or free energy at finite temperature) with respect to R. 

Another motivation is of a more physical nature, namely to investigate 
constraints for non-Newtonian gravity from the Casimir effect (cf., for in- 
stance, Refs. [28] and [29]). Present Casimir force experiments are so accu- 
rate that the possible influence from extra dimensions is taken seriously. The 
hypothetical extra force is usually taken to be of the Yukawa form. 

In the present paper we will however not consider possible Casimir- 
induced deviations from Newtonian mechanics, but focus instead on the anal- 
ysis of Casimir free energy and force in the presence of extra compactified 
dimensions at finite temperature T. At the final stages of our work, we ac- 
tually became aware of a series of papers of Teo, which go along similar lines 
as our considerations do. There is necessarily some overlap between Teo's 
papers and the present one. However, since our results were for the most 
part obtained independently of each other, it may be of interest to make a 
more close comparison: 

In Ref. [5], the finite temperature Casimir force in the presence of extra 
dimensions was calculated for the scalar field with Dirichlet boundary con- 
ditions. The free energy was not considered. The force was always found to 
be attractive. When the ratio R/L (in our notation) increased from 0.1 to 1 
the Casimir force was found to increase quite much, by a factor of about 3. 
In our approach, the Casimir force is calculated via the free energy. A sat- 
isfactory feature is that a detailed comparison (some calculation is needed) 
shows our results to be in agreement with each other. 

A more extended version of this work was given in the next paper [6], 
considering also the free energy. In the third paper in the series [7], the free 
energy was studied for a massive scalar field with general curvature coupling 
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subject to Dirichlet or Neumann boundary conditions, thus a more general 
setting than that considered in our present investigation. 

Finally, in three recent papers Teo investigated the finite temperature 
effect for a scalar field in multidimensional space assuming Robin boundary 
conditions [22]; then gave a comparison between Kaluza-Klein and Randall- 
Sundrum models with special emphasis on the non-reversal of sign in the 
presence of warped extra dimensions [23] , and worked out the Casimir theory 
for a multi-layer dielectric slab model in piston geometry |24j . 

In the present paper we focus on the following two issues: 

1) We allow there to be a 'refractive' index, called N, experienced by the 
scalar field in the physical medium between the plates. This influences the 
Casimir free energy as well as the entropy, and is of interest in connection 
with the Nernst theorem as T — > 0. To our knowledge, such a generalization 
has not been considered before. The introduction of a refractive index is not 
so trivial as one might imagine, if one compares with electrodynamic theory: 
in that case, the Casimir free energy expression F gets an extra factor N in 
the denominator at T = 0, the finite T corrections being more complicated. 
Here the T = behavior can be understood in a physical way by observing 
that the momentum of a photon in the medium can be written in the form 

+ n 2 m 2 /R 2 , where kj_ is the photon momentum parallel to the plates 
whose separation is R, and m is an integer (in turn, this is intimately linked 
to the assumption that the photon momentum is as given by the Minkowski 
energy-momentum tensor). Taking into account the volume element in phase 
space, one obtains the essential result for F. Cf. Ref. [23], and also the 
remarks of Ravndal [2SJ [27] . 

2) Our second point is to take into account the frequency dispersion 
of the medium in a crude way, by including the refractive index for low 
frequencies only. Specifically, we cut off the influence from the medium, thus 
setting N — 1, for all frequencies exceeding the first Matsubara frequency 
(i = 2ttT. From a basic viewpoint this procedure is permitted, since there is 
no conflict with the fact that the refractive index - like any other generalized 
susceptibility - has to be such that it decreases monotonically along the 
positive imaginary frequency axis. Moreover, in physical units £i is lying 
in the region 10 14 s _1 and is thus quite high, making the model sensible 
physically. This generalization is made for the case q = 1 only. 

We begin in the next section by considering the general case where the 
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number p of edges in physical space is allowed to take the values 1, 2, or 
3. Thereafter we specialize to the case p = 1, corresponding to the conven- 
tional setup with two parallel plates separated by a gap R. Zeta function 
regularization is employed. Our main focus will be on low temperatures as 
mentioned, as this is of main physical interest. For iV constant we evaluate 
the Casimir free energy, the corresponding pressure, the entropy, and work 
out some simplifying limiting expressions. The dispersion generalization is 
considered in section 6. 

Finally we mention for completeness that all plates, in physical space as 
well as in the extra compactified space, are so large that edge effects are 
negligible. We use natural units, with h = c = k B = 1. 

2 General formalism, when p = 1, 2, 3 

We assume finite temperatures from the outset. Let p be the number of edges 
in physical space. Thus p can take the values p = 1 (two plates), p = 2 (four 
plates), or p = 3 (box). The number M of transverse dimensions in physical 
space is M = 3 — p. We let index i± refer to the transverse directions, so 
that i± G [1,M]. The number q of extra compactified dimensions is at first 
assumed arbitrary. The total spacetime dimension D is thus D = 4 + q. 
We assume there to be a scalar field <&(x l ,t) in the bulk. As mentioned, 
Dirichlet boundary conditions are assumed on the plates in physical space. 
The plate separations are Ri, with i G In the extra space with di- 

mension q we assume, in accordance with common usage, that a torus of 
circumference 2irLj, j G [l,q], is attached to each spacetime point. Periodic 
boundary conditions are assumed for the extra dimensions. 

As mentioned, we will allow for a constant 'refractive index' N in the 
physical bulk. This index describes the response from the homogeneous and 
isotropic medium to the scalar field, in the same way as the usual refractive 
index yfe\x does in a dielectric medium, e being the permittivity and [i the 
permeability. Our model is thus a sort of scalar electrodynamics which has 
been considered occasionally in earlier studies also, for instance in Ref. [30J, 
dealing with curvilinear space. In the ordinary case D = 3 our model implies 
that the scalar photons possess the momentum k = Nu\t, in accordance with 
the Minkowski momentum in dielectric media in ordinary electrodynamics 
as already mentioned, implying a spacelike total photon four-momentum. 

The central differential operator in Euclidean space (r = it) is = 
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N 2 d 2 + J2i + dj. The free energy F can be found by using the zeta 
function for the operator — 

U(^) = EV, (i) 

J 

where in the present generalized case the eigenvalues are 



Here J is an index referring to all the indices {rii}, {mj}, {k t }, as well as the 
Matsubara frequencies (i = 2irlT with / = 0, 1, 2, .... The Dirichlet conditions 
on the plates causes the first sum to run over positive rii only. Further, rnj 
extends over all integers because of periodic boundary conditions. In fl2} we 
have also introduced k\ as the square of the transverse wave number k^, i.e., 
k\ = Yli ± =i • We emphasize that N refers to the physical space only. To 
introduce an analogous refractive index in the extra space would hardly have 
any physical meaning. 

The general expression for the free energy reads 

F = -^Co E (0). (3) 

We first consider the interior free energy, called Fj, in the cavity. With V± 
denoting the transverse volume, V± = Yli ± =i^i±i anc ^ with the quantity A 
defined as 

^« + £(f)\£(|)\ (4) 

i=l N ' j=l x J ' 

we have for the zeta function 

-r t OO OO OO « 

uw = (^mE E E jd M k x (ki+Ar. ( 5 ) 

(=— oo {rii}=l {m,j}=— oo 

Here the notation {rij} = 1 means that n\ = 1,2, ...oo, n 2 = 1,2, ...oo. The 
integral can be evaluated using the technique of generalized polar coordinate 
transformations [31]. This yields, after integrating over all angles, 



DC 



T _ M OO OO oo 

UW^f^E £ £ /. *r-V + A)~. (6) 

V 1 K 1 ' l=-oo {m}=l {mj}=-oo ' 
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After the variable change x = r 2 /A the integral part can be recognized as 
the integral representation of the beta function B{y,u) = T(v)T(u)/T(v + u) 
and hence the zeta function is 

l=— oo {m}=l {irij}=— oo 

The free energy can easily be found from the above expression if we remember 
that (g(z) /T(z))' z=0 = g(0) for any function g(z). 




Figure 1: Illustration of the four cavities of the piston model. 

We now introduce the piston in the i = p direction, observing that the 
Casimir free energy of the piston geometry, as usual in Casimir calculations, 
consists of four parts: Fi = Fj(Ri, . . . , R p -±, R p , L\, . . . L q ) is the free energy 
of cavity I with length R p in the piston direction. For cavity II we have 
Fji = Fu(Ri, . . . , Rp-i, X — R p , Li, . . . L q ) and similar expression for cavity 
III and IV with lengths X/rj and X(l — 1 / 77) in the piston direction respec- 
tively. The parameter 77 is taken to have a value of around 2 (see Fig. 1). 
The main point here is that cavity I is a squeezed region, in which the dis- 
cretization of allowable wave numbers is essential. Cavity II, on the other 
hand, is non-squeezed, so that summations over wave numbers can be re- 
placed by integrals. Similarly, the reference configuration (cavities III and 
IV) corresponding to rj ~ 2, is also non-squeezed. The technique is explained 
in more detail in Ref. [32J, Ch. 4, for instance. The conventional situation 
with parallel plates is retrieved by taking X — > 00. 

Thus in F n the sum over the discrete k p = n p 7i / R p goes to an integral in 
the limit X — > 00. This integral can be evaluated together with the integrals 
over the transverse directions in Eq. (jSJ). To find the free energy of cavity II 
we need only to replace M with M+l, V± with (X — R)V±, remove (n p Tr/R p ) 2 
in A and remove the sum over n„. 
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Assume hereafter that the compactified extra dimensions have all the 
same size, (Lj = L). We introduce the Epstein-like zeta function 



E v (a 1 ,a 2 , ...,a v ; s) = ^ \^2 a 

K}=i \i=i 



(8) 



After some algebra we then obtain the following expression for the free 
energy, called F, of the piston system 



F =F l + F n - F m - F w = --Tir M / 2 



(2tt 



fc=0 
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This expression holds for arbitrary temperature, for arbitrary p G [1,3] and 
for arbitrary integers q, when all the Lj are equal. As is usual, F refers to 
unit surface area. Recall that the dimensions of the plates in the transverse 
dimensions are assumed infinite. 



3 The case p = 1 

When p — 1, the system in physical space consists of two parallel plates 
separated by a gap R p = R. Thus M — 2. The value of q is at present kept 
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arbitrary. With T(-3/2) = (4/3) v^r we get 



T 

8^ 



E 

k=0 



(2) 



q—k 



rf-r 



7T 2 1 
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;-l 



+ 2 ((2vriVT) 2 , ^, i, . . . , i; -l) } 
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q—k 



1 1 

, ... , 
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(10) 



+ 2£ 1+g _ fc ((2vriVT) 2 , -L . . . , -L; -? J } 



The terms with k = q are independent of L and are equal to the Casimir 
free energy in ordinary 3+1 spacetime (q = 0). Naming these terms F q= Q we 
evaluate them separately 



9=0 



T 

8^ 



( ^ 



\R 2 ' 



r (-i) i ft I «; -i ) + 2«2 ( P^rAT)', w -, -1 



7T 



(11) 



Note that the term E q _ k . . . , -p-; — |j in Eq. (TTDjl vanishes when k = q. 
Equation (ITT]) can be simplified by using Ei(ai] s) = a]" s Ci?(2s), where (r is 
the Riemann zeta function, and the reflection formula 



Cr(z) = 



;i2) 



The Casimir free energy per surface area is then 

9-1 



T 

8^ 



E 
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+ 2£ 1+ ,_ fe ((2vriVT) 2 , -L . . . , -i; -~J } 



(13) 



with 



^=o = N 3 T^ - T^|- - T^±E 2 ( (2vriVT) 2 , --; -1 ) . (II) 



The Epstein-like functions in Eqs. ( TT3"j) and (|14|) will be zeta-function regu- 
larized [32j [33] . By repeatedly using 

E v (ai, a 2 , . ■ . , a v ; s) = —-E v _ 1 (a 2 , a 3) . . . , a v ; s) 

1 ^Trfs-i)^ / 1 
+ n\ T=^—\ — £V-i(a 2 ,a 3 , . . . ,a y ; s - -) 

2 y oi 1 (s) 2 

™ / t/ \ -fs-l/2)/2 



X ^ 1/2 ( ^k ni ( £ 



ni,ri2,...,nv=l \i=2 
V \ 1/2^ 



J=2 



we can express F in terms of the modified Bessel functions of the second 
kind, K v . We will now discuss limits of physical interest. 



4 Approximations 

As the above formalism is complicated, it is of interest to make mathematical 
simplifications applicable to situations of physical interest. 

There is one simplification that is obvious, resulting from the fact that 
the separations in physical space space are very much greater than those in 
the extra space. Casimir measurements are usually done with separations of 
order 1 /xm or less. In the extra space, one can estimate the typical distances 
by looking for minima of the free energy. For instance, by making for a 
moment a digression into warped space, in Ref. [MJ we found that at low 
temperatures there is a minimum of the free energy in a Randall- Sundrum 
model if the length r c of the orbifold is of order 10~ 18 GeV -1 ~ 10~ 19 fm. It 
appears reasonable that we can set 

R/L > 1 (16) 

for all practical purposes. 
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4.1 Low temperatures 

The most important limiting case is when T is low. This case gradually 
catches up the quantum mechanical zero-point energy as T — > 0. We shall 
take the low-temperature limit to mean 

RT < 1. (17) 

Obviously, it implies that also LT < 1. In dimensional terms, Eq. (II 7p 
implies (kBT)(R/hc) <C 1, or RT <C 0.002 with T in kelvin and R in meters. 
Thus to one percent accuracy, R = 100 nm corresponds to T < 200 K, 
whereas R = 1000 nm corresponds to T < 20 K. 

We now use Eq. (}T5]) to find the finite expression for F q=0 , 



q ~° 90 16R 2 8tt v ! 1 \R 2 

1 -eJ^-A 



12ttN 'V^R 2 ' 2 

3 



n,l=l ' 



2 D_2 



1440JVi2 3 90 v/2 V/? 3 ^ll/ *\NTRJ 



(18) 

Notice that ^T(— f^-; — lj cancels — T , coming from the Mat- 

subara integer equal to zero (I = 0), in view of Eq. (fl2"]l . This equation 
does so far not imply the low temperature limit. But since the argument 
in K v is large, we need only the first term of the asymptotic expansion 
K v {z) = (ir /2z) 1 / 2 e~ z . Due to the product nl in the exponent of the expan- 
sion we only need the dominant term n = I = 1 in the double sum. Thus 
F q =o is 

F q - = ^_ + N*T^ - iVT 2 ^ (19) 



in the low temperature limit. Here T 4 ^- is the Stefan-Boltzmann term cor 



90 

>2 



1U0NR 3 90 2R 

man 

responding to the vacuum energy in empty 3+1 spacetime and — e~N^r 
is the leading correction term. One may notice that there is no term propor- 
tional to T 3 £r(3); this term being characteristic for electromagnetic fields. 
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As regards the dependence upon N, the following points may be noted. 
At T = the (free) energy is one half the sum over eigenfrequencies u. In 
view of the Dirichlet boundary conditions on the plates this sum becomes 
transferable to a sum over integers, which may conveniently be done in view 
of zeta function regularization. As u = k/N, the proportionality factor 1/N 
in the first term in Eq. (fT9"j) necessarily follows. The second term in the same 
equation involves the thermal occupation number \j(e w l T — 1). This does 
not contain N in itself, but the refractive index occurs again in the phase 
space element d 3 k/(2n) 3 , which brings in a factor N 3 . This is seen to agree 
with Eq. (fT9|) . The other terms in Eq. (1TB]) are more complicated, but it is 
useful to observe that F is generally for finite T equal to T multiplied with 
the logarithm of the partition function. As iV occurs only in the combination 
(NT) in the quantity A defined in Eq. (j3J), which in turn is closely related 
to the eigenvalues of the basic Xj operator in Eq. (jSJ), the ratio F/T has to 
be a function of (NT). We see that this is the case for all terms in Eqs. ffT8l) 
and f|T9|) . Cf. also the remarks of Ravndal [261 127] . 

To find the low temperature limit of the part of the free energy depending 
on L and q it is convenient to evaluate the remaining part of Eq. (1T5|) in pairs. 
In this way we obtain after a lengthy algebra (details omitted) the following 
expression for the free energy in the low temperature limit 



T(-2)E g 




1 



+ N a T 



7T* 1 



L 2 



;-2 



(20) 



This is our main result. At zero temperature only the first three terms 
remain. The result agrees with Ref. [UJ when N — 1, except from terms in 
that reference leading to repulsive forces. We £1X6 ? clS mentioned, concerned 
with the piston model, where the repulsive terms are removed. This is also 
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the situation considered in Ref. [TU] . 

Again using Eq. ( |15p we can rewrite the zero-temperature energy E per 
unit area as 

2 <?-l , 



k=0 v 7 

J 1 1 1 3^ 



q—k 

2 



1/2N 



8L 2 Rtt 2 N ^ \^ 3 ) \ L , 

mi,m 2 ...,m,_ k ,n=l \i=l / \ \j=l 



(at T = the thermodynamic energy i£ is the same as the free energy F). The 
Epstein-like terms are independent of i? and will not influence the Casimir 
force. The arguments of the modified Bessel functions are proportional to 
R/L 3> 1 and will only give exponentially small corrections to the first term. 

The two last terms in Eq. ( 1201) are low temperature corrections from the 
compactified extra dimensions. They are proportional to the number q and 
are exponentially decreasing in 1/T. 

The Casimir force P per unit area can be found as P = —dF/dR. Since 
K' v {z) = -\[K v ^{z) + K v+1 (z)} we find, using Eq. (TU]), 



dF Q - 7T 2 , T ,^ 4 7T 2 1 m /iVT 

p g~ u N T T \ / 

9=0 dR A80NR 4 90 ^ U 3 



3n2 / 7rZn \ 7rn2 / / 7rZn \ / 7r/n \\ 

2^1 § ViVT^J " 2NTllF> V H^Ry U^JJ 



(22) 



2 2 

> N T h T p tttr 

480iV# 4 90 2F 3 



When iV = 1, this is one half of the result given in Eq. (3.16) in Ref. 
when dealing with the electromagnetic field. The factor of 1/2 difference is 
what we would expect, since we are considering the scalar field. The more 
general case with refractive index N was recently considered in Ref. [25J. 

We evaluate the other part of the free energy (pertaining to the extra 
dimensions) in the same way. Making use of Eq. ( 1T5|) for the Epstein-like 
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functions, we obtain the following final expression for the Casimir pressure 
(details omitted): 



71- 



9-1 

E 



480NR* 8L 2 R7T 2 N^\k 

k=0 



(2) 
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+ q R3 eXP \ NT\ R 2 ' L 2 



(23) 



When N — 1, this expression agrees with Teo [5] (she does not restrict all 
separations in the extra space to be equal). The T = terms arising from 
the compactified dimensions are negative, corresponding to an attractive 
force. The temperature corrections (the last three terms) are small in the 
low temperature limit. 



4.2 High temperatures 

We shall consider only briefly the case of high temperatures, 

RT > 1. 



(24) 



For instance, in order to satisfy this condition with one percent accuracy 
when the separation is R = 1000 nm, the temperature must be quite high, 
T > 200000 K. At such temperatures ordinary solid bodies do not exist. 
This limit is thus of less physical interest than that of low temperatures. 
Let us consider, however, the surface force density. Our starting point is 
again Eq. f fl3|) . from which we extract terms pairwise. Omitting all details 
we present only the final result for the high-temperature limit of the Casimir 
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pressure: 



p = _ rU3) _ 2nTW c _ 4nRNT qT 



SnR 3 R 2tiRL 2 ^ 

qT (1 + (2nNTL) 2 ) exp (^-2^ I + (2it NT L) 2 " 



tcRL 2 



When N — 1 our results are in accordance with Ref. [5] though here we have 
included exponentially small corrections as well, not just the terms linear 
in temperature. In Ref. [5] the size of each extra compacitifed dimension is 
assumed arbitrary. 



5 Entropy, and the Nernst theorem 

Of interest is here the low temperature limit. To evaluate the entropy S = 
—dF/dT we extract once again terms from the expression (Fl3|) pairwise. As 
above we abstain from giving the details of the calculation and present only 
the result: 
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(26) 



(terms that do not contribute in F are omitted). As the derivatives of K v 
will produce terms containing K v with the same argument we see that, when 
T — > 0, the K v will decay exponentially implying that S — > 0. The Nernst 
theorem is satisfied. This is as we should expect, in the present case with 
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idealized boundary conditions. The same property is known to hold for a 
metal without extra dimensions, when ideal boundary conditions are assumed 
for all frequencies. (The current discussion about thermal Casimir corrections 
relates to the case where dissipation is included; for a discussion on these 
issues see, for instance, Ref. [36].) 

6 On dispersive properties, when p = q = 1 

As the second novel element in our analysis we will give a simplified analysis 
of the change in free energy due to frequency dispersion in the medium filling 
the physical space. As before, we assume p = 1 so that there are two infinite 
plates in physical space separated by a gap R. To simplify the formalism 
we assume moreover that q — 1, so that there is only one extra dimension 
(containing a vacuum, not a medium). This formal simplification will still 
allow us to show the characteristics of the dispersive behavior. 

Now write the free energy F as a sum of its q = part (i.e., ordinary 
space, with no extra dimension) and a q = 1 part, 



Here the first term is given by Eq. (ITT]) , for arbitrary temperature. The 
Matsubara frequencies are d = 2ttIT with I = 0, 1, 2, ... The first finite value, 
corresponding to I = 1, is in dimensional terms Ci — ^ksT/h = 2.4 x 
10 14 s -1 . This value is quite high, and it is therefore physically reasonable, 
as mentioned in the Introduction, to ignore the dispersive properties for 
higher values of I. We accordingly adopt in the following a dispersive model 
in which N(() = N = constant for I — 1 in the Epstein-like function E 2 , 
and then falls off abruptly such that we may put N(Q = 1 for / > 2. This 
point simplifies the formalism, but still allows us to see the characteristics of 
the dispersive behavior. 

Starting from the general expression for E 2 we can thus make the following 
effective substitution 



F 



F q=0 + F i 



q=l- 



(27) 
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where 

00/ v \~ s 

E^( ai ,a 2l ...,a v ;s) = ^ c2 + E a ^i) ( 29 ) 

{n 3 }=\ \ j=l J 

is a generalized Epstein-like zeta function. The first term in Eq. ffTTj) {F q=0 ) 
does not refer to the dielectric properties of the physical space. As for the 
first term and all other terms not containing N(() we get no change from 
the previous sections. 

Then, for the third term with E\ in Eq. (|TT|) we make use of the same 
effective substitution as above, extracting the I = 1 term, 

Ex ((2vriVT) 2 ; -£) E, (W) 2 ; -£) + (2vrT) 3 (iV 3 - 1) (30) 

No regularization of the iVo term is here necessary. Combining all three terms 
in F q=0 we get 



-<n=i T ( „2ttTN f^ 2 . \ P 2xT ( 71-2 . 



+ -^(2vrT) 3 (iV 3 -l). 



(31) 



Here we have collected all terms that give the free energy without disper- 
sion, i.e. is equal to Eq. ( ITT1) with iV = 1 for all I. Hence, we can 
use all expressions obtained above and only concentrate here on the addi- 
tional terms introduced by the dispersion model. Since E 2nTNo (7r 2 /R 2 ; —1) = 
n 2 / R 2 (eh((2RTN ) 2 ; — 1) where (eh is the Epstein-Hurwitz zeta function, 
we might recall the conventional expansion 



(eh(p;s) = ^2(n 2 +p)' 



71=1 



V* 1 v^r(s-i/2) /2 

2 P 2T(s) y 



+ — ^ £ nS ~ 1/2 Ks-x^n^p)- (32) 

^ ' 71=1 

cf. [33J (p. 81). Although this series is for many purposes a convenient 
expression, it is in our case with small arguments better to write (eh i n 



-s/2+1/4 



oc 
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another way: as — s is an integer we make use of Eq. (3.24) on p. 37 in 
Ref. E2. 1 to write 



E c 1 (l;s) = ( EH (c 2 ;s) 

K _ LJ A k + a )^ mu , o^Jk ^ r(„ ^ j-a. (33) 

k=0 

With this we find the finite expression 
Ttt 



; T(-l)( EH ((2RTN y;-l] 



AR 2 ' 

= ~ % (r(-l)C^(-2) - (-£) (2RTN )^ (34) 

= - t -^At*nI 

8nR 2 3 

We can now insert these expressions into Eq. (J3il to evaluate F q=0 at ar- 
bitrary temperature. Moreover, to simplify the formalism somewhat, we 
restrict ourselves to the low temperature approximation, RT <C 1. From 
Eq. (USD, d3U) and (JMD we find 



F q - = — RT\N*-l)-^- + ^-T> e -^, RT«1, (35) 
9-0 3 v o ; 1440.R 3 90 2R ' V ; 

where the influence from the medium is shown explicitly in the first term. 

Consider next the term coming from the extra dimension, q — 1. From 
Eq. (TI3D we first get 

RT „ ( \ 3\ T „, _ „ /tt 2 1 



RT / 1 3\ T / 7r 2 1 

+ _ B2 ( (arJW) . )--r(-i)^((ar«T)» (-10 



3 Equation (3.24) in [32] p. 37 gives three different expressions for the Epstein-Hurwitz 
formula. One is for — s ^ N and 1/2 — s ^ N, the second for — s £ N, and the third 
for 1/2 — s € N. According to our recalculations two of these expressions are correct, 
but in the case (— s £ N) the expression should be different and as given above. See the 
Appendix. It should also be noted that an assumption for these expressions is that c 2 is 
small. This assumption is fulfilled in the low temperature limit since c 2 = (2RTNq) 2 <C 1. 
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where the terms are written in order of increasing complexity. The two last 
terms depend on N and can be rewritten in the same manner as Eq. (I28p . 
Again, we write the free energy as the free energy of a medium with N = 1 
(for all I) plus correction terms from 1 = 1, 



-H-l) ( E^ Na ( i; -l\ - E? T { 4; -1 



(37) 

From p. 48 in [33J we adopt a formula similar to Eq. (fl5|) (with some re- 
writing) 

Ey(di, a 2 , • • • , a v ; s) = — -E v _ 1 (a 2 , a 3 , . . . , a v ; s) 

i nrr(s-l) i 
+ ^- r( s ) ^~ l(a2 ' fl3 ' ' • • ' av; s ~ 2) 

27r * -0+l/2)/2 \^ s-1/2 / 2 , 2^\ ( 38 ) 



x K s _i /2 [ ^=nt ( c 2 + ]T 




i=2 



Then, we can write 
T( i\F 2wTNo (— —■ 1 



+ 



m,', 



O 00 / 2 \ 3/4 / I 2 

_ ^ n _ 3/2 ,m_ + {27tTNq) 2\ Ks/2 2 Rn^ + (2nTN ) 2 

^,n=l ^ \ 



(39) 

There occurs a cancellation since -T/(2n)(R/(2 v ^))T(-3/2)E 2wTN °(l/L 2 ; -3/2) 
-RT/(3n)E 2wTNo (l/L 2 ; -3/2). With the aid of Eq. ([33]) and the reflection 
formula for (r we see that 

l-TV 1 \ Tp2nTN ( 1 i | Cr(3) 8 4 3 3 



-_r(-l)^— _!j =-^ + -vr 4 T^. (40) 
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When all is put together, we get 
4 
3 



F q=1 = F q N = f - -7r 3 T 4 L(iV 3 - 1) 



T + (27rTN )A ^ K 3/2 ( 2RnJ ^ + {2irTN Q ) 2 



^ U 2 V V ' V V L 2 

m,n=l x 



+ ^= f: + (2.T) 2 ) ^ K 3/2 I 2^ ^ + (2*1? 

m,n=l 



(41) 



Since R/L > 1 we only need to take into account the n — m — 1 term of 
the double-sum. Some calculation yields for the medium-independent F^ 1 - 
terms 

r(-l)^a ((2.T) 2 , ± -l) = -\V { -1)E 2 -L -1 

1_ 3gCfi(5) 1 y> m 2 /2nmiA 

720TL 3 32vr 5 TL 4 2ixRTL 2 ^ n 2 2 \ L ) 

m,n=l x ' 



2V2^T , m/^ 2 " 2 «i 2 \ /' /A 2 in , 

+— E l rVir + i?l- (42) 



and 



1 3\ 1 9 



2>T ^ m 2 //m\ , inS 



i,m=l 



At last we have the final low temperature expression for L/i? 1 



T 2 /tt 2 1 1/4 
W 1 <= tVj? 



n \ R 2 + L 2& 



2nR 

T 
2~kR 



+ (2nTN ) 2 e- 2R ^ +im)2 
^ + (27tT) 2 e- 2R ^ +(27TT)2 , RT < 1, L/R < 1. 



(44) 
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Adding Eq. (1351) and (|44p we get the total free energy per unit surface in 
physical space to the lowest order in RT and R/L 



1440ttL 3 



~R 



90 V 7T 



L)-7r 3 T 4 (iV 3 -i) 



T 



,e at 



T 2 /tt 2 



e " 2R / L — —a/- h — 

2i? 8(nRL) 3 / 2 n \ R 2 L 2 



(45) 



T 



2vri? V L 2 



^ + (2vrTiV ) 2 e- 2 V^ + (™) 2 



+ 



2vri? V L 2 



4 + (27rT) 2 e~ 2i ViM 2 ' rT)2 , #T < 1, L/i? < 1. 



For the pressure P = —dF/dR we obtain, when omitting the small correction 
terms, 



P 



TT 



2 n 2 T 4 



+ T 



480i? 4 90 

— — i 

e at 1 



[240 (iV 3 - 1) + 1] 



-2H/L , 
+ 



1 M,l 



(46) 



2i? 3 AL{nRLf/ 2 

- — (\ + (2nTN )A e ~ 2 ^^TN 0) > 
ttR \L 2 ) 

+ ? ( J_ + (2VTT) 2 ) e "W£^ 

Here the first term gives the conventional attractive force in physical space. 
The second term is attractive and can be traced back to the q = con- 
tribution. In the force expression there are Bessel functions with the same 
arguments as in the free energy. There are several sums over Bessel func- 
tions with different arguments; we have kept the largest term from each sum 
to show the relative sign and size. Note that since iVo > 1 the contribution 
from the extra dimension is attractive. It is seen that cross terms, containing 
contributions from both the R and the L spaces, are all exponentially small. 



7 Summary 

We began by developing the finite temperature Casimir theory of the scalar 
field, assuming Dirichlet boundary conditions for two large parallel plates 
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with a gap R in physical space endowed with a medium with constant 're- 
fractive index' N, and in addition an arbitrary number q of extra compact- 
ified dimensions. Zeta function regularization of the Epstein-like functions 
was used throughout. In the main part of the formalism the radii L of the 
compactified dimensions were taken to be equal. For general values of T, 
the free energy F was given by Eq. ( TT5j) . For low T (i.e., RT <C 1) which 
was the case of main physical interest, the approximate expression for F 
was given by Eq. (120]) . and the corresponding Casimir pressure was given by 
Eq. ( 123|) . When N — 1, agreement was found with Teo [5]. The Nernst the- 
orem was found to be satisfied; cf. the low temperature entropy expression 
(12 6p . The case of high temperatures was briefly discussed. In all cases of 
physical interest, the inequality R/L ^> 1 could be assumed. 

In the dispersive part of the theory considered in Section 6 (assuming 
q — 1), we introduced a rough model implying neglect of the dispersive prop- 
erty at all frequencies higher than the first Matsubara frequency, £i = 2nT. 
It is worth noticing that for the Epstein-Hurwitz zeta function we found it 
convenient to use the expansion (155|) , instead of the more conventional expan- 
sion (152~j) . We found (155|) to be the most suitable form for low temperatures. 
For RT < 1 we calculated the free energy in Eq. (H5|) . showing contributions 
from physical space as well as from the extra dimension. The correspond- 
ing pressure in Eq. ( 146|) contains both attractive and repulsive terms. The 
common property of the cross terms, i.e. terms containing both R and L 
contributions, is that these terms are exponentially small. 
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A Regularization of the Epstein-Hurwitz zeta 
function 



In this appendix we will argue why the correct regularized expression for the 
Epstein-Hurwitz zeta function should be 



E c 1 (1;s) = Ceh(c 2 ;s) 



-s-l 



k=0 



k\T(s) 



( R (2k + 2s)c 



2k 
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1/2 



2T(s) 



1 

s 

2 



A-2s 



(A.l) 



for — s G N. As noted above, this expression is a bit different from the one 
given in Eq. (3.24) in [32] : 



E^s) 



E(-l) fc£ i^C*(2fc + 2 S )c 2fc 



k=0 



klT(s) 



(A.2) 



It is seen that the two expressions differ with respect to two terms: the first 
is k = — s from the sum over k, the second is a term proportional to c 1_2s . 
Referring to the basic calculation in [32] we will argue why the k = —s term 
should be omitted and the c 1_2s term should be present. 
We start with (3.7) in [32] 



E1(l;s) 



k V ( k + S J, „ , 2A 



k=Q 

kjtl/2-s if 1/2-seN 
k^-s if -seN 



k\T(s) 



( R (2k + 2s)c 2k + Res 2 , (A.3) 



where 
Res, 



Res 



T(s + a) 

r(a + i)r( s ; 



c a Ci?(2s + 2a)7T esc 7ra; a 



(A.4) 



for z = 1/2 — s and z = — s. (The function 0(a) comes from the asymptotic 
behaviour of the integrand when the sum is converted into an integral; cf. 
[32].) The sum over k runs from to oo, but skips k = 1/2 — s (if 1/2 — s G N) 
to avoid the pole Cr(1) an d k = —s (if — s G N) to avoid the pole T(0). The 
next step is to calculate the residue at a = 1/2 — s and a = —s. We quote 
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the results from [32] : 



Resi/ 



2-s 



r (s - 1) 



■2s 



2r( s ) J 

r( fl )r(§- a ) 

"1 f , fl\ , /3 



- I 2~ S) +1 ° g C 1+7 



where ip(a) = 4- logT(a). For a = — s we have 



Res_ 



if 1/2 - s £ N, 



if 1/2 -seN, 



(A.5) 



|c~ 2s if - s i N 
if - s G N. 



(A.6) 



The results can be divided into three: 1/2 — s N and — s N; 1/2 — s G N; 
and — s G N. We are interested in — s G N. Combining Eqs. ( IA.3j) . ( 1A.5|) and 
( ]A.6j) we find this to be equal to Eq. (lA.ip . The sum over k only runs to 
— s — 1 since all terms with k > —s are equal to zero since Cr( — 2n) = for 
n = 1,2,3... and k = — s is excluded for — s G N. The c 1_s term must be 
included since — s G N implies 1/2 — s ^ N. 
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